Chiral Magnetic Effect and Chiral Phase Transition 
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We study the influence of the chiral phase transition on the chiral magnetic effect. The azimuthal 
charge-particle correlations as functions of the temperature are calculated. It is found that there is a 
pronounced cusp in the correlations as the temperature reaches its critical value for the QCD phase 
transition. It is predicted that there will be a drastic suppression of the charge-particle correlations 
as the collision energy in RHIC decreases to below a critical value. We show then the azimuthal 
charge-particle correlations can be the signal to identify the occurrence of the QCD phase transitions 
in RHIC energy scan experiments. 
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The phase transitions of quantum chromodynamics 
(QCD), for example the evolution between chiral sym- 
metry breaking and its restoration, the color dcconfine- 
ment and confinement, have been one of the most active 
topic in nuclear and particle physics in recent years [lj. 
Such phase transitions can be driven by the temperature 
and density of the system. It is then expected that these 
phase transitions occur and the deconfined quark gluon 
phase (QGP) is formed in ultrarelativistic heavy-ion col- 
lisions 01 (for example the current experiments at the 
Relativistic Heavy Ion Collider (RHIC) and the upcom- 
ing experiments at the Large Hadron Collider (LHC)) 
and in the interior of neutron stars [4-6]. However, the 
explicit variation behavior of the signals to identify the 
phase transitions with respect to the temperature and 
density needs further investigations. 

Recently, The STAR Collaboration at RHIC report 
their measurements of azimuthal charged-particle corre- 
lations in Au + Au and Cu + Cu collisions at ^/snn = 
200 GeV. They find a significant signal consistent with 
the charge separation of quarks along the system's or- 
bital angular momentum axis 0, H[ . The observed charge 
separation indicates that parity-odd domains, where the 
parity (V) symmetry is locally violated, might be cre- 
ated during the relativistic heavy- ion collisions [il-fllj]. 
The charge separation is related with the so called "chi- 
ral magnetic effect" which means that a magnetic field 
in the presence of imbalanced chirality induces a current 
along the magnetic field, and therefore results in that 
positive charge is separated from negative charge along 
the magnetic field [111 ]. 

Now that the chiral magnetic effect can be observed 
through the measurements of azimuthal charged-particle 
correlations in the relativistic heavy-ion collisions, a nat- 
ural question arises, i.e. whether can we detect the prop- 
erties of the QCD phase transitions, especially the chiral 
phase transition through the observations of the chiral 
magnetic effect? To answer this question, we have to 



study how the chiral magnetic effect or the charge sepa- 
ration effect is influenced by the chiral phase transition. 
This is our central subject in this letter. 

In this work, we will study the chiral magnetic ef- 
fect and the QCD phase transitions in the 2+1 flavor 
Polyakov-Nambu-Jona-Lasinio (PNJL) model [12|. The 
validity of the PNJL model has been confirmed in a series 
of works by confronting the PNJL results with the lattice 
QCD data [IM!]. The PNJL model not only has the 



chiral symmetry and the dynamical breaking mechanism 
of this symmetry, which are same as the conventional 
Nambu-Jona-Lasinio model, but also include the effect 
of color confinement through the Polyakov loop. There- 
fore, the PNJL model is very appropriate to describe the 
QCD phase transitions at finite temperature and/or den- 
sity. 

The Lagrangian density for the 2+1 flavor PNJL model 
is given as 
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D n = d n _ iAI ± with A n = s ^ A o^ A a _ gA° a ±f = 
—iA^. X a are the Gell-Mann matrices in color space; 
to = diag(mQ , Wq, mg) is the three-flavor current quark 
mass matrix. In this work, we take TOq = TOq = to , while 
keep TOq being larger than m l . U(<&, $*,T) in the PNJL 
Lagrangian density is the Polyakov-loop effective poten- 
tial, which is expressed in terms of the traced Polyakov- 
loop $ = (Tr c L)/N c and its conjugate = (Tr c Zt)/jV c . 
In this work, we use the Polyakov-loop effective potential 
which is a polynomial in $ and $ 

M($,$*,T)_ &a(T) 



13(, given by 



T 4 



-$* 3 )+^($*$)2( 2 ) 



2 



with 



62 (T) = a + ai 



«2 



«-3 



(3) 



Parameters in the effective potentiai are fixed by fit- 
ting the thermodynamical behavior of the pure-gauge 
QCD obtained from the lattice simulations. Their val- 
ues are ao = 6.75, a\ = —1.95, 02 = 2.625, 03 = —7.44, 
63 = 0.75 and 64 = 7.5. The parameter To is the 
critical temperature for the deconfinement phase tran- 
sition to take place in the pure-gauge QCD and To is 
chosen to be 270 MeV according to the lattice calcula- 
tions. Furthermore, we also need to determine the five 
parameters in the quark sector of the model, which are 
771 q = 5.5 MeV, m s = 140.7 MeV, GA 2 = 1.835, KA 5 = 
12.36 and A = 602.3 MeV. They are fixed by fitting 
ttt-tt = 135.0 MeV, m K = 497.7 MeV, m v > = 957.8 MeV 
and f n — 92.4 MeV [lB|- 

In the parity-odd domains which are created dur- 
ing relativistic heavy-ion collisions, the number of left- 
and right-hand quarks is different because of the axial 
anomaly. In this work we introduce the chiral chemical 
potential 115 to study the left-right asymmetry following 
the method of Ref. [Hj], where the chiral chemical po- 
tential /U5 is related with the effective theta angle of the 
9- vacuum through ^5 = do8/2Nf and Nf is the num- 
ber of flavor. Consequently, we should add the following 
term 
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to the Lagrangian density in Eq. (p}, where /15 = 
diag(/^g, /xf, Next, we consider the case that a ho- 
mogenous magnetic field B is along the direction of the 
orbital angular momentum of the system produced in a 
non-central heavy-ion collision. In the following we de- 
note this direction with z-direction and particle momen- 
tum in this direction with p%. 

In the mean field approximation, the thermodynamical 
potential density for the 2+1 flavor quark system under 
a homogeneous background magnetic field B and with 
left-right asymmetry is given by 
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where 



e/l = \^^\q f \eB+p 2 , 
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with qi(i — u,d,s) being the electric charge in unit of 
elementary charge e for the quark of flavor i and the 
constituent mass Mj reading 
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and (4>ip)i is the chiral condensate. In Eq. (J5J) we also in- 
clude the quark chemical potential The momenta of 
charged particles in the longitudinal direction, i.e., the z- 
direction, are not influenced by the background magnetic 
field and ps in the expression of the thermodynamical 
potential density in Eq. (|S|) is continuous; while the mo- 
menta in the transverse plane are discretized due to the 
magnetic field effect. |e/| in Eq. ([6]) is similar to the mag- 
nitude of the momentum in free space, s (for fermion and 
for anti-fermion is — s) in Eq. ([5]) is the helicity of particle 
and we should emphasize that at the lowest order of the 
transverse quantum number, i.e., 77 = 0, the quark spin 
only has one value in the z-direction, which means that 
charged particles in the lowest transverse level are polar- 
ized by the external magnetic field; however particles in 
higher levels, i.e., 77 > 0, are not polarized. Therefore, 
the charge separation effect only comes from quarks in 
the lowest transverse level. 

In order to relate our calculations with observable in 
heavy-ion collisions, we define A + (A_) to be the positive 
(negative) charge difference in unit of e (— e) between on 
each side of the z = plane, which is also the reaction 
plane. Here we use the notations in Ref [ltj. Taking 
particles with positive elementary electric charge e for 
example, we can express A + as 

A + = J d 3 a#7°7V)n=o 

= J d^^Vii - ^7Vi)n=0, (9) 
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The subscript n = in Eq. ([§]) indicates that only the 
lowest transverse level states contribute to the A + . In the 
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same way, we can obtain A± for the 2+1 quark system. 
We take A + for example once more, which is given as 
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eB 



dp?,— 



f(E u 



d P3TT f(E d + (l d - TTM, 



-f(Ed 



E d 

+ W/ 5 > 



E d ' 



dp: 



P3_ 

'E s 



f(E s 



P3 



E s 



(11) 



where V is the volume of the system, Ef is given by 
Eq. with n = 0, and 
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In fact, we can also obtain Eq. (fTTj) through differen- 
tiating the thermodynamical potential in Eq. ([5]) with 
respect to the chiral chemical potential (i\ and sum- 
ming over the contributions from positive quarks or anti- 
quarks. 

In the same way, differentiating the thermodynamical 
potential with respect to the quark chemical potential 
Hi and summing over contributions from the three-flavor 
positive quarks or anti-quarks we obtain the total positive 
electric charge number N + in unit of e, i.e., 
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Similarly, the total negative electric charge number iV_ 
in unit of — e can also be obtained. 

In the experiments of heavy-ion collisions, the az- 
imuthal charged-particle correlations, i.e., (cos(</> Q + 
4>p — 2^ bp)), are used to detect the ^-violating ef- 
fect 0, S, Efj. Here </> and ^ rp are the azimuthal an- 
gles of the particles and reaction plane, respectively, a, 
j3 represent electric charge + or — . With the notation 
a a/3 = —{cos((f) a + (f)p—2 i ^jip)}, it can be shown that [1(3] 
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where the azimuthal angle distribution of the charged 
particles is assumed to be 



dN± 1 „ 1 . . 
do 2n 4 



(17) 



Since we mainly focus on the influence of the QCD 
phase transitions, especially the chiral phase transition, 
on the chiral magnetic effect in this work, we will ne- 
glect the screening suppression effect due to the final 
state interactions [1 01 ] and make fa = 0, then we have 

a ++ = a = — a -| Therefore, we only study a ++ in 

the following. 

Minimizing the thermodynamical potential in Eq. ([5]) 
with respect to three-flavor quark condensates, and 
$*, we obtain a set of equations of motion. We neglect 
the influence of the magnetic field on these equations of 
motion in our numerical calculations, since the magnetic 
field (eB = 10 2 ~ 10 4 MeV 2 in the non-central heavy- 
ion collisions [l(|) has little impact on these equations of 
motion. 
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FIG. 1: Correlation a++ as function of the temperature calcu- 
lated in the PNJL model with (15 = 150 MeV (left panel) and 
(is = 250 MeV (right panel). The magnetic field corresponds 
to eB = 5 x 10 3 , 10 4 , and 5 x 10 4 MeV 2 from top to bottom, 
respectively. 



In Fig. [T] we show a ++ defined in Eq. (TPo)) as function 
of the temperature at several values of the chiral chemical 
potential (15 (here (15 = (i§ = (i§ = and the magnetic 
field strength. We find that there is a pronounced cusp in 
a ++ at the critical temperature during the chiral phase 
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transition (the critical temperature T c = 209 MeV for 
fi s = 150 MeV and T c = 185 MeV for /i 5 = 250 MeV in 
the PNJL model). From the Fig. [T] one can also find that 
although the value of a ++ is proportional to the square 
of the magnetic field strength, the shape of the curve for 
a ++ as function of temperature is almost independent of 
the magnetic field strength. Furthermore, the cusp at the 
critical temperature in the curve becomes much sharper 
with the increase of the chiral chemical potential. With 
the decrease of the temperature, when the temperature 
is below T c , chiral symmetry is dynamically broken and 
quarks obtain masses. Since the axial anomaly can be 
suppressed by the mass effect, which has been discussed 
in detail in Ref . 17J j the chiral magnetic effect can also be 
suppressed by large constituent quark masses. Therefore, 
the azimuthal charged-particle correlations described by 

a ++ (a and ) defined in Eqs. (|15[) (|16p are quite 

decreased once the temperature is below the critical tem- 
perature. It can been seen from Fig. [I] that, when the 
temperature is above T c , a++ decreases with the increase 
of the temperature, which is because higher temperature 
makes it more difficult to polarize quarks with magnetic 
field and thus suppresses the charge separation effect. 

What do our calculated results imply in future energy 
scanning experiments of heavy-ion collisions? With the 
decrease of the heavy-ion collision energy, the tempera- 
ture of the QGP produced in the fireball at early stage 
is also decreased. Since the magnetic field produced in 
non-central collisions decays rapidly with time (Io| . the 
observed charge separation mainly carries the informa- 
tion of the QGP at early stage. Therefore, we expect 
that the azimuthal charged-particle correlations (espe- 
cially for the same charge correlations, because the op- 
posite charge correlations are suppressed by final state 
interactions) increase as the collision energy is lowered. 
However, when the collision energy is lowered to the 
value that cannot drive the chiral phase transition, it 
is expected that the azimuthal charged-particle correla- 
tions are quite suppressed. Therefore, we can employ the 
charge separation effect to locate where the QCD phase 
transitions occur. 

In summary, we have studied the influence of the QCD 
phase transitions on the chiral magnetic effect. The az- 
imuthal charge-particle correlations as functions of the 
temperature are calculated in the PNJL model. It is 
found that there is a pronounced cusp in the azimuthal 
charge-particle correlations around the critical tempera- 
ture of the chiral phase transition. We predict that there 
will be a sudden suppression of the charge-particle corre- 
lations with the decrease of the collision energy. It indi- 



cates that azimuthal charge-particle correlations can be 
a signal to identify chiral phase transition in the energy 
scan experiment in RHIC. 
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